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1 INTRODUCTION 



ABSTRACT 

In this paper we describe the Bayesian link between the cosmological mass function 
and the distribution of times at which isolated halos of a given mass exist. By as- 
suming that clumps of dark matter undergo monotonic growth on the time-scales of 
interest, this distribution of times is also the distribution of 'creation' times of the 
halos. This monotonic growth is an inevitable aspect of gravitational instability. The 
spherical top-hat collapse model is used to estimate the rate at which clumps of dark 
matter collapse. This gives the prior for the creation time given no information about 
halo mass. Applying Bayes' theorem then allows any mass function to be converted 
into a distribution of times at which halos of a given mass are created. This general 
result covers both Gaussian and non-Gaussian models. We also demonstrate how the 
mass function and the creation time distribution can be combined to give a joint den- 
sity function, and discuss the relation between the time distribution of major merger 
events and the formula calculated. Finally, we determine the creation time of halos 
within three N-body simulations, and compare the link between the mass function 
and creation rate with the analytic theory. 
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as the epoch at which non-linear collapse is predicted. Two 
derivations were given, one o f which directly used the trajec 



> The hi erarchical build-up of self-gravitating dark matter is tories invo k ed in PS theory dPeacock & Heavens 1990 



though to drive evolution in the observable universe. The 
formation ol clumps ol dark matter precipitates the lorma- 
tion of galaxies by providing a potential well into which 
gas can fall and subsequently cool. Violent mergers between 
equally sized halos and their associated galaxies are thought 
to be important for starbursts and quasar activation. In or- 
der to model and understand the observable universe it is 
therefore essential to understand the build-up of the dark 
structure. 

The most widely used analytic model for the distribu- 
tion of mass in isolated halos at any epoch come s from Press- 
Schechter (PS) theory ( [Press fc Schechter 1974 ). By smooth- 
ing the initial field of density fluctuations on different scales, 
information on the distribution of perturbation sizes can be 
obtained. Linking the time at which these perturbations col- 
lapse to the initial overdensities using the simplified spher- 
ical top-hat collapse model allows the distribution of mass 
in isolated halo s at any epoch to be de t ermined ([Press & 



Bond 



Schechter 1974 



Peacock fc Heavens 1990; Bond et al. 1991) 



In Percival & Miller ( 199E ) (hereafter paper I), we used 
the tenets of PS theory to model the related, but distinct 
problem of determining the distribution of times at which 
halos of a given mass are created. Here, 'creation' is defined 



et al. 1991 ), and one of which used Bayes' theorem to con- 
vert from the PS mass function to a time distribution. The 
second derivation required the prior for the creation time 
which was calculated by examining the trajectories model. 

In this paper we extend the Bayesian link between the 
mass function and the creation time distribution to cover 
any mass function. This is important, not only because it is 
known that standard PS theory is wrong in detail (e.g. Sheth 
& Tormen 1999), but especially because the new extension 
applies to mass functions derived from more general density 
fields including non-Gaussian models (e.g. Matarrese, Verde 
& Jimenez 2000). 

First, we adopt the assumption that all clumps mono- 
tonically increase in mass on the cosmological time scales of 
interest. This monotonic growth is an inevitable aspect of 
gravitational instability. Every epoch should now be thought 
of as a creation time for a given clump, and we need not 
make the distinction between the creation time distribution 
and the distribution of times at which a given halo exists. 

In order to convert from a mass function to a distribu- 
tion in time we require the prior for the creation time. This 
is the rate at which creation events occur, given no informa- 
tion about the halo mass. In this work we use the spherical 
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top-hat collapse (STHC) model to provide a simple mech- 
anism for determining the required rate. In Section ^ we 
derive the link between collapse time and the overdensity at 
an early epoch for the STHC model within any Friedmann 
cosmology. Having determined that this relation is indepen- 
dent of halo mass, this leads directly to an approximation 
to the prior for the creation time, described in Section |j| 
This is the second major assumption adopted in this paper: 
that the prior for the creation time is well approximated 
by this simple model for the break-away of structure from 
linear expansion. This means that following the two simple 
assumptions detailed above, we are able to convert any mass 
function to give the distribution of epochs at which halos of 
a given mass are created. 

Simple models of cosmologically evolving phenomena 
often adopt an important mass range rather than a specific 
halo mass (e.g. paper I, Granato et al. 1999). In order to use 
the work presented here in these models, the joint distribu- 
tion of halos in mass and creation time is required. Although 
calculating the required joint probability is formally impos- 
sible because the equations cannot be properly normalised, 
a formula with the correct shape can be determined and is 
presented in Section |j| 

So far we have not made a distinction between the slow 
accretion of mass onto a halo and major mergers between 
halos. Such a distinction is important because only ma- 
jor mergers are thought to play a vital role in starbursts 
and quasar activation (see paper I). The time distribution 
calculated in this paper determines when halos existed (or 
were created by any mechanism assuming monotonic clump 
growth) which is not necessarily equal to the distribution of 
merger events. This is discussed in Section ^. 

Finally, we compare the analytic link between the mass 
function and the creation rate to the results from three nu- 
merical simulations of structure formation in different cos- 
mological models. An analytic fit to the mass function as 
described by Sheth & Tormen ( 1999 ) is adopted and is con- 
verted into a creation rate using the STHC model. This 
model is compared with and shown to be in good agreement 
with the numerical results. 



2 THE SPHERICAL TOP-HAT COLLAPSE 
MODEL 

In this Section we analyse the STHC model which is the 
simplest model for the way in which clumps of dark mat- 
ter break free from linear growth and undergo non-linear 
collapse. We present the derivation of the link between the 
initial overdensity and the collapse time t co u in a form which 
clearly shows that this link is independent of the mass of the 
overdensity. The derivation also demonstrates a method for 
calculating this link within any Friedmann cosmology. Sim- 
ilar derivations have been previously discussed in a variety 
of subsets of this space: for an Einstein-de Sitter model, 
a derivation is given by Gunn 
£lv = model by Lacey & Cole i 
universe by Eke Cole & Frenk 
results is given in Kitayama 




Gott (1972), for an open 
and for a flat flv 7^ 
A sum mary of these 



Suto ( 1996 ). A numerical 
prescription for the calculation of t he overdensity in any 
cosm ology has also been developed (Somerville & Primack 
1993). 



These derivations all use the same basic idea which is 
adopted in this work: the behaviour of two spheres of equal 
mass is compared within the cosmological framework. One 
of the spheres evolves with the background density pb{i), 
while the other is perturbed by a uniform excess density 
Ap(t). In subsequent analysis, a subscript denotes that a 
quantity relates to the sphere with background density, and 
'p' to the perturbation. 

Matter is assumed to be an ideal fluid with no pres- 
sure and the universe is modelled as spherically symmetric 
around the perturbation. As a consequence of Birkhoff 's the- 
orem, the gravitational field of both the perturbation and 
the background is described by a Robertson- Walker (RW) 
metric with curvature constant K, and RW scale factor a(t). 
The behaviour of such perturbations is governed by Fried- 
mann's equation which we will consider in the form: 



'da\ 2 2GM . 2 v 2 



(1) 



where M is the mass inside the sphere. Note that in order 
to compare spheres with different behaviour, we do not nor- 
malise the scale factor a(t) to equal the curvature scale (by 
dividing by so K is allowed to take any real value. 

To calculate the behaviour of the overdensity at an early 
time, we note that a series solution for a(t) in the limit 
t — > can be obtained for Equation |l|. This is given by 
a = at 2 ' 3 + + 0(t 6 ? 3 ), where: 



_ /9GA/\ 1/3 



(2) 



' 20 \GM) 

Using the fact that the spheres contain equal mass, the be- 
haviour of 8(t) = Ap(t) I pb{t) in the limit t — » is given 
by: 



Mm (Wf) 

t-o \ p b (t) J 



a p (ty* 



= -(/3 p -/3 6 )^ /3 + 0(t 4/3 ). 
a 



Denning: 



K 



(3) 



(4) 



(GM# ) 2/3 ' 

the present day normalisation of Equation [j] gives that for 
a sphere of uniform background density: 

/ 2 \ 2 / 3 

e b = (n M + Clv-l) (ts— , (5) 

We can now combine Equations ^, |^ & |E] to determine the 
behaviour of 8(t) in the limit t — ► as a function of e p : 

9 /4\ 1//3 
Urn S(t) = — (-) 
t-a v ; 20 \3J 



+ O[(H t) 4/:i ] 



2/3 



(Hot) 



2/3 



(6) 



If the field of perturbations is linearly extrapolated to 
present day and nor malise d here, the approximation of Car- 
roll, Press & Turner ( 1992 ) for the ratio of the current linear 
amplitude to the Einstein-de Sitter model can be used to 
extrapolate the limiting behaviour of S to this epoch. The 
extrapolated limit, (5ii m , is related to e p by: 
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Shr, 



\2/3 



(o M + n v - i) 



4/7 



fiv + 1 



2 
1 



2/3 



1 + ^v 



(7) 



A similar formula is possible if the field of fluctuations is nor- 
malised at any other epoch. Note that 5u m tx (e p + constant) 
and the time dependence of <5n m is given by that of e p . 

For the perturbation, the radius of maximum expansion 
can be calculated from Equation |l|: this radius corresponds 
to the first positive root of the equation 2GM + HqQ.vo?' — 
Ka — denoted by a m ax- This leads to a necessary and suf- 
ficient condition for the perturbation to collapse: that such a 
(finite) root exists. Because of the symmetry in Equation |l| 
this model predicts that the perturbation will collapse to 
a singularity at a time equal to twice the time required to 
reach maximal expansion: 



Hot c 



-1/2 



da* 



(8) 



where we have changed from a to a* = aH^ 3 /(GM) 1 ^ 3 , 
and aJn ax is the first positive root of the equation: 



2 + fl v {a*) - e p a = 0. 



(9) 



Although collapse to a singularity does not occur in practice, 
the virialisation epoch is assumed to be similar to i C oii- 

For perturbations that collapse, <5u m is called the 'criti- 
cal' density and is denoted S c - Equation ^ then gives 5 c (e p ). 
Equations]^ & [] give t C oll(ej>), and the combination of these 
three Equations gives the required link between <5 C and the 
collapse time. Note that these Equations are independent of 
the perturbation mass, and therefore so is the link between 
the initial overdensity and the collapse time. 

In practice we wish to use these Equations to calculate 
8c(z co n,QM,flv) or d8c{z co \i,Q,M,Q.v)/dt where z co \\ is the 
collapse redshift. Unfortunately this is not easy as Equa- 
tions ^ & |^ cannot be inverted to give e p (t C oii). The proce- 
dure adopted is as follows: the collapse time can be numer- 
ically determined from z co \\ using the Friedmann equation 
for the background cosmology. e p can be determined nu- 
merically using Equations ^ & [], and S c can be calculated 
using Equation [?]. dS c /dt can be calculated numerically from 
^c(tcoii) and is discussed further in the next Section. 

For the subset of cosmological models with Q.v = 0, the 
above procedure is simplified and analytic formula can be 
obtained for 8 C . In this case, Equation M reduces to: 



ri/'p , 2 . 
Hoicoli = 2 J y — — e p J 



-1/2 



da* 



(10) 



Making the substitution tan(#) = (2/a* — e p ) 1//2 , this inte- 
gral can be solved to give: 

2tt 



Hot c 



3/2 ' 



(11) 



We now sho w tha t these Equations provide the result 
of Gunn & Gott (1972) for an Einstein-de Sitter cosmology. 
In this case, substituting Equation |ll] into Equation |^ gives 
that: 



(12) 



We can now change from collapse time to collapse redshift 
to give: 



8c(z coll ) = ^(127r) 2/3 (l + ^coii) 
~ 1.69(1 + z con ), 



(13) 



which is the equation of Gunn & Gott (1972) 



3 FROM A MASS FUNCTION TO A TIME 
DISTRIBUTION 

In this Section we show how to convert from a mass function 
to the distribution of times at which isolated halos of a given 
mass exist. First, we make the assumption that the mass of 
any clump is a monotonically increasing function of time so 
that the mass will increase between any two epochs. This is 
true for Press-Schechter theory (see paper I) . Note that this 
mass growth is not constrained to be continuous and the 
mass is allowed to undergo instantaneous finite increases, 
or 'mass jumps'. Following this assumption, every epoch at 
which a halo exists should also be considered as a 'creation' 
epoch: every halo is a new isolated halo of some mass. The 
distribution of 'creation events' is therefore the same as the 
distribution of times at which the halos exist. Note that by 
definition this only applies to isolated halos which have not 
been subsumed into larger objects. 

In this paper we have called this epoch the 'creation' 
time of a halo in order to avoid confusion with other au- 
thors definitions of the 'formation' time of a halo. Note that 
this semantic change was not adopted in paper I. Th e 'for - 
mation' time of a halo was defined by Lacey & Cole (199;) 



as the latest time when the largest progenitor of a halo has 
a mass less than half that of the final halo. This definition 
makes sense if we are discussing a non-evolving quantity, say 
the existence of a galaxy halo, and wish to know when it was 
formed given that it exists at present day. However, suppose 
we do not know anything about the build-up of a halo be- 
fore or after it has mass M and only wish to know when 
it was likely to have existed. This Lacey & Cole definition 
of 'formation' cannot help us for we do not know the time 
and mass from which to determine progenitors: progenitors 
of what? 

In order to calculate the probability density function 
(pdf) of the times at which halos exist, we consider the set 
of all possible times and all possible halo masses. This is the 
'sample space' of our 'experiment'. The experiment consists 
of choosing a particle, or small mass element, and an 'event' 
is given by any subset of the sample space: for instance that 
the particle is part of a halo of mass Mi < M < M2 created 
at time t\ < t < t%, or that the particle inhabits a halo of 
mass M, created at time t. 

Denoting a generic pdf by the function /, the mass func- 
tion is given by f(M\t) dM, the distribution of halo masses 
at a given epoch. This is equal to Mn(M)/p where n(M) is 
the number density of halos. The pdf we wish to determine 
is given by f(t\M) dt, the distribution of times at which ha- 
los of mass M were created. Note that our assumption of 
monotonic mass growth means that t is the same variable 
in f{M\t) dM and f(t\M) dt. These pdfs are then related by 
the following formula, based on Bayes' theorem: 
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f(t\M) dt = 



f(M\t) dMf(t)dt 
f °°[f(M\t)dMf(t)}dt' 



(14) 



where f(t) dt is the normalised prior for time, or the distri- 
bution of creation events in time given no information about 
the mass of halo. In paper I we calculated the prior using 
the Brownian random walks invoked in PS theory with a 
sharp fc-space filter. In order not to bias the distribution of 
up-crossings within this model, we assumed a uniform prior 
for 5 C . 

The reason the prior is uniform in S c follows from the 
STHC model. Within this model, it is the density S associ- 
ated with a particle that is important, and the barrier has 
to move from S to 5 — dS for 'halo creation' to have occurred. 
Given that the mass of all clumps monotonically increases, 
all particles will be associated with creation events at any S c . 
Following these two observations, any two equal width in- 
tervals in 5 C should contain equal 'numbers' of halo creation 
events. 

The derivation presented in the previous Section showed 
that for the STHC model, the link between the critical over- 
density and the collapse time is independent of the perturba- 
tion mass. Therefore, given no information about the mass 
contained within a perturbation, the pdf for the time at 
which the perturbation collapses should be assumed to be 
proportional to the time derivative of S c (t). This gives the 
rate at which the collapse threshold 8 c {t) crosses the initial 
overdensities. This can be calculated numerically from the 
following formula: 



dS c 



de v 



dt X dt 



d 
de v 



( a max( e p) 



— + S7y(a*) 2 — e p ) tin 



V ° 



, (15) 



where a^ ax is the first positive root of the equation 2 + 
Qv(a*) 3 — e p a* — 0. Note that for cosmologies with Q,v = 0, 
the above equation can be analytically solved as for Equa- 
tion [ll], and the derivative dS/dt is proportional to t~ 5 ^ 3 . 

Unfortunately, dS c (t)/dt cannot be normalised so that 
it integrates over all time to give unity. This means that we 
cannot simply take a multiple of dS c (t)/dt as the prior for 
the collapse time. However, we can still use Equation [l4| by 
making use of a mathematical trick and placing an arbitrary 
upper limit on t, t u , which can be removed later without 
affecting the result. This gives that: 



f(t\M)dt= lim 



f(M\t) dM f(t,t u )dt 
jjf" [f(M\t)dMf(t,t v )] dt 



(16) 



The connection between the mass function and the cre- 
ation rate of halos presented in this Section is consistent 
with that of paper I: the prior for time used is exactly the 
same. We have merely shown that adopting the STHC model 
for the rate at which structures are created allows any mass 
function to be converted to give the pdf of the time at which 
a halo of a given mass is created. 



(17) 



4 THE RELATION WITH THE 
MULTIPLICITY FUNCTION 

Changing variables from mass to a function of 

v= — 

om 

alters the form of the standard PS mass function to one 
which is invariant with respect to time. Here cm is the rms 
fluctuation of the initial density field smoothed with a top- 
hat filter on a scale related to mass M. Unless stated oth- 
erwise we change variables in the mass function from M to 
\nv(M,t). The normalised pdf f(\nv\t) is called the multi- 
plicity function and is related to the mass function by: 

d In v 



f{M\t)=Af{\nu\t) 



dM 



(18) 



where A is a normalisation constant. Note that we have re- 
tained the condition on time in /(ln^|t), to emphasise that 
we are still concerned with the distribution of halos at a 
particular epoch. Although the multiplicity function has a 
form which is invariant with respect to time, it still gives 
the distribution of \xiu we would expect for halos given a 
particular time. This is not the same as the the distribution 
of In v we would obtain if we chose halos at random in both 
mass and time, or the distribution of \av wc would obtain 
if we chose halos only of a particular mass. 

If the mass function can be written in a form which is 
independent of time as described above, then under the same 
change of variables, the creation rate becomes independent 
of halo mass. The resulting pdf /(lnf|M) is now only valid 
if we are examining the distribution of halos at fixed mass. 
Following the notation adopted above, this is given by: 

d In v 



f(t\M)=Af(lnu\M) 



dt 



(19) 



5 THE JOINT DISTRIBUTION OF HALOS IN 
MASS AND TIME 

Because the mass of each halo is assumed to monotonically 
increase with time, within any interval of mass and time, an 
infinite number of 'creation events' occur. This means that 
the joint probability of the existence of a halo in both mass 
and time cannot be properly normalised. 

Equation |l^ gives the link between two pdfs, the mass 
function and the creation rate using a mathematical trick 
to cope with an un-normalised prior in time. The numer- 
ator of this equation is the joint distribution of halos in 
mass and time, f(M, t) dM dt = f(M\t) dM f(t) dt. The de- 
nominator is not a function of time: it only normalises the 
resulting formula so f(t\M) dt integrates to unity. Follow- 
ing this argument, given a mass function, multiplying by 
d5 c (t)/dt creates a function with both the correct mass and 
time behaviour. This joint distribution function (not a pdf) 
has the same mass dependence as f(M\t) and the same time 
dependence as f(t\M). 

As an e xamp le we consider the fitting function of Sheth 

6 Tormen (1999) to the multiplicity function determined 
from the results of N-body simulations for different cosmo- 
logical parameters: 



Mn(M) 



dM = f(M\t) dM = /(In u\t) din v 
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V^( i+ i)^' 2/2dini/ ' 



(20) 



where and a & p are parameters. Note that Sheth 

& Tormen displayed this formula using a different notation 
to that adopted here, although parameters a and p are the 
same in both cases. A is determined by requiring that the 
integral of /(In v\t) over all In v gives unity. Sheth & Tormen 
found best fit parameters a = 0.707 and p — 0.3 for their 
simulations and group finding algorithm. The standard PS 
multiplicity function has a = 1, p = and A = 1/2. Un- 
less stated otherwise, by standard PS theory, we refer to the 
adoption of this multiplicity function combined with top-hat 
filtering (to calculate erf./). In order to convert this function 
to provide a model of both the time and mass of halo cre- 
ation events, all we need to do is to multiply by dS c /dt. 

For standard PS theory, writing v explicitly in terms of 
<j%[ and 5 C we find that the joint distribution of the existence 
of a halo in mass and time reduces to: 



f(M,t) dMdt = 



(27T) 



x exp 



2< 



da* 



dM 



dS c 
~~dt 



dM dt, 



(21) 



Although not normalised, such a formula integrated over any 
two areas of the mass-time plane will provide the correct 
relative number densities. 

Note that this is not the same formula as obtained by 
simply multiplying the mass function with the creation time 
distribution at fixed mass. This would be inconsistent within 
a Bayesian framework and would produce a joint density 
function which lacks the correct mass and time behaviour: 
the form of each conditional pdf is altered by the other. 
Care should therefore be taken when using the creation rate 
in models which also include the mass function. 



6 THE RELATION WITH MERGER EVENTS 

So far, we have only been concerned with the epoch at which 
a halo is created. However, there is an important distinc- 
tion between major mergers and the slow accretion of mass 
when applying the results in models of certain cosmologi- 
cal phenomena. For instance, only violent merger events are 
thought to be important for starbursts and quasar activa- 
tion. In paper I, we showed that for standard PS theory with 
a sharp fc-space filter, if mass jumps in a particular trajec- 
tory correspond to merger events, then the distribution of 
mergers is the same as that of the build-up of matter from 
all types of creation event. This is because the trajectories 
are Brownian random walks which have the special property 
that their form is independent of the initial point. 

Given only the mass function and the assumptions out- 
lined above it is not possible to determine how each clump 
increases in mass, only the distribution of times at which it 
reaches a certain mass. More information about the build- 
up of individual clumps is required before the distribution 
of major mergers can be determined. Such information is 
available in PS theory and follows from the argument that 
each trajectory gives the history of the halo masses in which 
a particular small mass element resides. 
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Table 1. Table showing the parameters of the different cosmo- 
logical models adopted in the three N-body simulations. 



7 DESCRIPTION OF THE NUMERICAL 
SIMULATIONS 

A direct approach to modelling structure formation is to 
simulate the evolution of the mass density of the Universe 
using a distribution of softened particles. We have run three 
such s imulations using the Hydra N-body , hydrodynamics 
code (Couchman, Thomas & Pearce 1995) with 128 3 dark 
matter particles to model the build-up of halos for three 
different cosmological models, described in Table |l|. 

In order to determine the rate at which halos are created 
within these simulations, we output particle positions at a 
large number of times. For the PCDM simulation, we output 
particle positions at 362 different epochs, separated by ap- 
proximately equal intervals in time. For the OCDM simula- 
tion the number of outputs was 345 and for the ACDM simu- 
lation, 499. The box size chosen was 100 /i _1 Mpc for all three 
simulations which gave a particle mass of 2.6 x 10 11 M Q for 
TCDM and 7.9 x 10 10 M Q for the other two simulations. 
Groups of particles were found for each output using a stan- 
dard friends-of-friends algorithm with linking length set to 
b = 0.2 times the mean interparticle separation. 



8 FITTING TO THE MASS FUNCTION 

The multiplicity function averaged over all output times is 
presented from each of the simulations in Fig. [l]. Here we 
have only considered groups containing over 45 particles in 
order to limit the number of false detections due to numer- 
ical effects. In compiling the data in this way, we have as- 
sumed that converting from mass to In v does indeed convert 
the form of the mass function into one which is independent 
of epoch. This Figure has been produced in such a way as 
t o be directly comparable with figure 2 of Sheth & Tormen 
( 1909| ). For comparison we also plot their best fit model and 
the predictions of standard PS theory. 

W e have also plotted the model of Sheth & Tormen 
( 1999| ) (Equation |2o| ) after allowing the parameters to vary 
to simultaneously fit the data from all three simulations. We 
find slightly different best fit parameters to those of Sheth 
& Tormen. Our best fit parameters are a = 0.774, p = 0.274, 
compared to standard PS theory a = l,p = and Sheth & 
Tormen a = 0.707, p — 0.3. Note that the difference between 
our best fit values and those of Sheth & Tormen could be 
explained by the different group finding algorithms used. 



9 COMPARISON BETWEEN THE ANALYTIC 
AND NUMERICAL HALO CREATION 
RATES 

Although we have argued that the monotonic increase in 
mass means that all epochs are 'creation' times for a given 
halo, we cannot simply compare the creation rate formulae 
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Figure 1. The measured multiplicity function from each of the 
three simulations analysed for all groups containing over 45 par- 
ticles (solid triangles plotted with Poisson error bars). See Sec- 
tion ^ for details of these simulations. In compiling these data in 
this way we have assumed that the transfer of variables from M 
to ln^ where v = Sc/ctm does make the distribution of_masses 
independent of epoch as found by Sheth & Tormen (1999). For 
comparison we have also plotted the best fit model to their data 
(dashed line) and the predictions of standard PS theory (dotted 
line). The solid line shows the best fit m odel t o our data allowing 
the parameters of the Sheth & Tormen ( 1999| ) fitting function to 
vary. 



with the distribution of halo numbers at different epochs: 
each halo should only be counted once. To determine the 
distribution of creation times of halos of mass M, we there- 
fore sequentially analysed the FOF output from z — 50 to 
present day. All halos of mass > M were examined at each 
epoch to determine whether they were 'new'. The definition 
of 'new' adopted was that at least half of the particles in 
a halo were not included in any halo of mass > M at a 
previous output time. The number of these halos in the re- 
quired mass range was taken to be the minimum number 
which could have been created between that output time 
and the previous one. In order not to miss creation events 
where a halo was created and subsumed into a larger halo 
all within the time interval between two outputs, we anal- 



ysed the progenitors of all new halos with mass greater than 
the required range. Those with a progenitor distribution at 
the previous step which could sum to a halo of the required 
mass were recorded as a possible halo of the required mass. 
In this way we determined the minimum and maximum mass 
which could have been created in each time interval between 
output from the simulation. 

In Fig. |^ we plot the creation rate for halos within two 
narrow mass ranges. In order to obtain the maximum num- 
ber of creation events, we have used relative low numbers 
of particles in each group. Data are plotted for groups of 
between 45 — 50 and 100— 110 particles. These distributions 
are compared with the three multiplicity functions plotted 
in Fig. converted into creation rates by multiplying by 
dS/dt for halos of mass equivalent to 45 or 100 particles. 
These curves have been normalised to the low redshift data. 

All of the models reproduce the decrease in creation 
events to present day seen in the simulations. As output 
from the simulation occured after approximately equal in- 
tervals of time, the high redshift data suffers as the intervals 
contain relatively more creation events. This means that we 
cannot precisely follow the build-up of the clumps, and the 
difference between the maximum and minimum mass which 
could have been created in each bin is increased. This is par- 
ticularly noticable in the OCDM simulation where halos are 
created at earlier times and we have fewer outputs from the 
simulation. 

However, there is evidence that the solid line (calculated 
from the best fit to the mass function) also fits the creation 
rate data the best out of the three models plotted. As a 
rough guide to this, the root mean square value between 
the plotted data points and the model is 3.4 for this curve, 
compared to 6.7 for standard PS the ory, and 5.0 for the best 
fit model of Sheth & Tormen fll99S| ). Note that the form of 
the creation rate is strongly dependent on the parameter a 
in Equation and only weakly dependent on parameter p. 
This is consistent with the importance of these parameters 
for the mass function: parameter o controls the position of 
the high-mass cut-off, whereas parameter p controls the low- 
mass tail of the distribution. 



10 CONCLUSIONS 

We have demonstrated a simple method for linking any mass 
function to the corresponding distribution of times at which 
isolated halos of a given mass are created. In order to provide 
this link we adopted the assumption that the time scales of 
interest are those over which the mass of every clump can 
be thought of as monotonically increasing. The prior for the 
collapse time was estimated using the STHC model which 
ties in directly with PS theory, although the method does 
not use any of PS theory beyond that of the STHC model. 
We have presented a new derivation of the link between the 
collapse time and initial overdensity for this model which ex- 
plicitly shows that this link is independent of the halo mass 
and is applicable in any Friedmann cosmology. Multiplying 
the mass function by a function with no mass dependence 
and proportional to the time derivative of the critical over- 
density then provides a joint density function with the cor- 
rect behaviour for the creation of a halo in mass and time. 
Integrating over the resulting joint density function will give 
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Figure 2. The creation rate of halos as determined from the three cosmological simulations (solid circles) for a) 45-50 particles in a halo 
and b) 100-110 particles. These group sizes correspond to masses of aproximately a) 1.2 X 10 13 M e for TCDM and 3.6 X 10 12 M e for 
the two other cosmologies, and b) 2.6 X 10 13 Mq for TCDM and 7.9 X 10 12 Mq for the two other cosmologies. The error bars show 
the binned minimum and maximum mass which could have been created in halos of the required masses in each time interval. Symbols 
are plotted half way between the two. Note that the errors bars do not include counting errors. For comparison we also plot the three 
models of the mass function multiplied by d8 c /dt with lines as in Fig. [l| The models have been normalised to the low redshift data. 



the correct relative number densities of halos within different 
mass and time intervals. 

We have extended the analysis of N-body simulation 
results presented in paper I to cover three simulations of 
the build-up of dark matter within different cosmological 
models. Rather than using PS theory, we have demonstrated 
how a fit to the mass function may be converted to give a 

crpatinn ratp Out, of the thrpp fnnrtinns wp have compared 



to the mass function data, the best fit model for these data 



when converted to a creation rate also fits the creation rate 
data the best. This gives us confidence that the formalism 
presented here is sound, and should give accurate results in 
more general situations, m particular non-Gaussian models. 
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